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ABSTRACT 



The ultra-compact Low Mass X-ray Binary (LMXB) X1916-053, composed of 
a neutron star and a semi-degenerated white dwarf, exhibits periodic X-ray dips 
with variable width and depth. We have developed new methods to parameterize 
the dip to systematically study its variations. This helps to further understand 
binary and accretion disk behaviors. The RXTE 1998 observations clearly show 
a 4.87d periodic variation of the dip width. This is probably due to the nodal 
precession of the accretion disk, although there are no significant sidebands in the 
spectrum from the epo ch folding search. From the negative superhump model 



(ILarwood et. al. Ill996l ). the mass ratio can be estimated as g = 0.045. Combined 



with more than 24 years of historical data, we found an orbital period derivative 
of Porb/Porb = (1.62 ±0.48) X 10~^yr~^ and established a quadratic ephemeris for 
the X-ray dips. The period derivative seems inconsistent wit h the prediction of 
the st andard model of binary orbital evolution proposed bv iRappaport et. al. 
(119871 ). On the other hand, the radiation-driven model (ITavani et. al. Ill99ll ) 



may properly interpret the period derivative even though the large mass outflow 
predicted by this model has never been observed in this system. With the best 
ephemeris, we obtained that the standard deviation of primary dips are smaller 
than that of secondary dips. This means that the primary dips are more stable 
than the secondary dips. Thus, we conclude that the primary dips of X1916- 
053 occur from the b ulge at the rim instead of the ring of the disk proposed by 



Frank et. al. I (119871 ) 



Subject headings: stars: individual (X1916-053) — X-rays: binaries — X-rays: 
individual (X1916-053) 
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Introduction 



Some X-ray binaries show periodic dips in their X-ray hght curves. To date , there are 
eleven X-ray binaries that known for periodic dipping features (iRitter &: Kolb I l2003l ) . It 
is widely believed that these X-ray dips are caused by X-rays from the region around the 
compact object being periodically absorbed by the vertical structure on the outer part of 
accretion disk. The periodic dips provide us an opportunity to measure the orbital period 
and investigate the accretion disk dynamics of the binary system. 

The low-mass X-ray binary X1916-053 (4U1915-05), composed of a neutron star and a 
white dwarf companion, shows various timing properties. For t he binary orbita l perio d, it 



White fc Swank 



can be obtained by per iodic X-ray intensity dips in light curves. I Walter et. al. I (Il982l ) and 



White fc Swank 



( 1982 ) first discovered its X-ray dips with Einstein and OSO-8 observations. 
(Il982l ) found two possible periods of 49.93 ± O.OGmin and 5 0.06 ± 0.03min, 
which is the shortest one among all X-ray dippers. iHomer et. al. I (120011 ) used the 0-C 
(Observed - Calculated) techn ique on the 1998 RX TE observation data to calculate the best 
fitting period of 3000.6 ± 0.2s. Ichou et. al.~l J200lh used 1979-1996 data to derive the orbital 
period of 3000.6508 ± 0.0009s and established an X-ray dip ephemeris. In addition to the 
~3000s recurrent dips, secondary dips are also occasionally seen at phases from ~ 0.4 to 
~ 0.6 relative to the center of primary dips. 

The optical modulation period of X1916-053 is 3027.5510 =fc 0.00 52 (s), which is only ~ 



1% significantly longer than the X-ra y dip period (IChou et. al. Il200lh . Optical observations 



have revealed the 3000s X-ray period (iGrindlay Ill989l . ll992l ). whereas th e 3028s optical perio d 



with a series of ~3.9d side bands was detected in X-ray light curves (IChou et. al. Il200ll ). 



The optical period has ones been considered as the orbital period of the system and the X- 
ray dip period is the beat period of bin ary orbital period and a th i rd co mpanion with orbital 



period ~3.9d (IGrindlay et. al. Ill988l ). However, IChou et. al. I (120011 ) concluded that the 



X-ray dip period, like other dipping sources, is the orbital period. Further, they concluded 
that the optical modulation is likely caused by the orbital motion coupled with a 3.9d disk 
apsidal precess i on pe riod as the superhump in SU UMa type dwarf novae. Furthermore, 
Retter et. al. I (120021 ) discovered a negative superhump with a period of 2979.3s , which is 
the beat period of orbital period and ~4.8d disk nodal precession period. 



On the other hand, dip shape modulation was first found in Ginga data (lYoshida et. al 



19951 ). Due to s hort observati o n dur ation, the modulation period could only be constrained 
to ~ 5-6 days. IChou et. al. I (120011 ) found that the dip phase probably modulated within 
a 4.85d or 3.76d period by analyzing 10-day consecutive RXTE observation in 1996. A dip 
phase variation with amplit ude ~0.05 and period of 4.74d was seen in 10-day consecutive 



RXTE observations in 1998 ( Homer et. al. 



20011 ). 
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Dip shape modulation provides an opportunity to study the dynamics of the disk. Un- 
fortunately, dip shape has never been quantitatively investigat ed before. For examp le, in 



order to investigate the possible periods between 3.5 and 5d, iHomer et. al. I (120011 ) have 
considered the area of dip, length of dip, depth of dip, and phase of deepest dip, but none 
of them have clear indication of ~4d cycle. In this paper, we have developed a new method 
to parameterize the properties of the X-ray dip to further investigate its variation using 
archived data from various observatories with a total time span of more than 24 years. 

In section [2], we briefly introduce the historical data applied in this study and the 
corresponding basic reductions. Section [3] describes data analysis methods, including the 
definition of the dip center, width and strength, as well as the results, such as the periodic 
variation of dip width, the long-term orbital ephemeris and the stabilities of primary and 
secondary dips. Finally, we discuss interpretations from our analysis of the nodal precession 
of accretion disk and the period derivative of the system in Section HI 



2. Observation and Data Reduction 

Many X-ray observatories have observed X19 16-053 since its discovery. More than 24 
years of observation data, from OSO-8 (1978) to XMM- Newton (2002), have been reduced 
through standard processes and archived on the HEASARC (NASA's High Energy Astro- 
physics Science Archive Research Center) website. Table [1] lists the historical observations 
of X1916-053 used in this paper. RXTE performed two 10-day consecutive observations in 
1996 and 1998, which enable us to study dip parameter modulations over a time scale of 
several days. Simultaneous observations by RXTE and XMM-Newton were also proceeded 
in May 2002. 

X1916-053 was observed by XMM-Newton on 2002 September 25 for about 15000 sec- 
onds. The XMM-Newton data were reduced through the standard process with SAS (XMM- 
Newton Science Analysis System). All the instruments, including MOS detector, pn detector, 
RGS spectrometer, and OM module, were used for this observation. Although pn detector 
provided larger count rate and better signal to noise ratio, it only detected four complete 
dips while MOS detector observed five dips. In order to collect more number of dips for bet- 
ter statistics, we used MOS data for dip phase analysis. The barycenter correction, which 
corrects the time system from satellite to the barycenter of the solar system, is important for 
the timing analysis of such a short orbital period. SAS provided the barycenter correction 
subpackage, barycen, for the XMM-Newton data. After applying the barycenter correction, 
the data can then be used for further dip analysis. 
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X1916-053 was observed by ASCA on 1993 May 2 for about 67ks and by ROSAT on 
1992 October 17 for about 18ks. Since FTOOLS provided the tools for correcting the time 
system, we downloaded the event files for ASCA and ROSAT. The ASCA data barycenter 
correction was produced by the subpackage timeconv while the ROSAT data used bct+abc. 
In addition, the event files of BeppoSAX observations were included in our analysis. The 
observation time was approximately 82ks in 1997 and 98ks in 2001. The data were barycenter 
corrected with the FTOOLS subpackage earth2sun, which can only correct the effect of earth 
to the solar system barycenter. 

The data collected by RXTE, Ginga, EXOSAT, Einstein, and OSO-8 were downloaded 
from HEASARC in light curve format already reduced through corresponding standard pro- 
cesses. Additional column BARYTIME in RXTE light curve data contains barycenter cor- 
rected time system by FTOOLS subpackage fxbary. All the other data sets were barycentric 
corrected by earth2sun only. 



3. Data Analysis 



3.1. Dip Parameters 

Researc hers have long been aware of variations of dip phase and width , as well as other 
parameters (lYoshida et. al. Ill995l : IChou et. al. Il200ll : iHomer et. al. Il200ll ). Unfortunately, 
dip profiles are usually complicated, and finding a unique model of paran aeters applicable 



to all dips is difficult. Taking the definition of a dip center as an example, lYoshida et. al. 



(I1995I ) defined dip center time as the midpoint of selected "start" and "end" times of a dip, 
but the valu e highly depends on how to choose the "start" and "end" times of a dip. On the 
other hand, IChou et. al. I (I2OOII ) defined dip cei iter time by fit t ing a dip with a quadratic 
curve around the point with minimum intensity. iHomer et. al. I (120011 ) fitted the dip with a 
Gaussian profile. However, none of these methods are suitable for dips with complex profiles, 
such as Fig. [H It is therefore necessary to extract dip properties with a new method that is 
not only independent of dip boundary selection but which can also be applied to all kinds 
of dips profiles, regardless of their complexity. 

An X-ray dip is caused by the absorption of the bulge in the accretion disk. We can first 
divide a light curve into dip states and persistent (non-dip) states by roughly guessing the 
boundaries of dips (see Fig{T]). The observed count rate during the persistent state Iq when 
the bulge is far from the line of sight is estimated by fitting a straight line of neighboring 
persistent states. Much like calculating the gravitational center of an object, dip center time 
(or phase in the folded light curve) is the average time of the dip weighted by the difference 
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of the dip state and the predicted persistent count rates: 



TV 



5^(/o - li) ■ ti 



1=1 



N 



1=1 



where the ti and Jj are the time and dip count rate of i*^ bin, respectively. To evaluate the 
physical width of the dip, the dispersion is defined as: 



W 



\ 



N 



^(/o - li) ■ {U - tc 



i=l 



N 



(2) 



i=l 



which is equal to the second moment of the dip state. 

Analogous to the equivalent width in spectrographic analysis, the equivalent width of the 
dip, which can also be considered the strength of dip, is 



(3) 



1=1 



By the definitions above, as long as the boundaries reside in the persistent states neigh- 
boring the dip, the points which lay on the persistent state beside the dip make little con- 
tribution to the dip parameters (Jq — /i ~ 0). The parameters are therefore insensitive to 
boundary selection for calculations. Furthermore, these three values are well-defined regard- 
less of the dip profile complexity. This method can also be applied to the eclipse source. The 
only constraint is that a completely observable dip or eclipse is required. In order to check 
if the parameters are sensitive to the bo undary choice, a tes t was performed to a dip which 



is ha rd to fit with a quadratic function (jChou et. al. Il200ll ) or a Gaussian (IHomer et. al 
20011 ) as shown in Figure El We set the left boundary (i=l, see Figure [1]) at a fixed value 



and gradually change the right boundary (i=N) from dip state to persistent state. The test 
result showed that the derived parameters are insensitive to the right boundary as long it 
well resides the persistent state. Similar result was obtained with a fixed right boundary 
at persistent state and a variable left boundary. We therefore conclude that our method is 
adoptable to dips of various profiles. 
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Simultaneous XMM- Newton and RXTE observations provide an opportunity to inspect 
whether the dip parameters are energy (or instrument) dependent. To test energy depen- 
dency, the XMM-Newton events detected by pn detector were divided into to a soft band 
(0.7-1.7 keV) and a hard band (1.7-8.0 keV) by the counts available roughly equally. Tabled 
lists the parameters of four dips in pn observation. There is no significant energy dependency 
in dip dispersion, but the equivalent width is highly energy dependent as expected since the 
dip is caused by absorption. The dip center time seems to have little energy dependency. 
This test was also performed with RXTE PCA data (2-9 keV) and XMM-Newton simulta- 
neous observations. Dispersions measured from both instruments are almost the same, but 
the energy dependency of dip center times is still present (see Table [3]). The difference is 
apparently due to the asymmetric bulge structure with different optical depths in different 
energy bands. Although dip center times are energy depend ent, the differences a re only ~ 



0.01 cycle, much smaller than the phase jitter (~ 0.05 cycle, IChou et. al. I (1200 ll )) and can 
be neglected in long-term analysis. 



3.2. Periodicity of Dip Parameters 

The new dip parameter definitions can be utilized to systematically study the variations 
of the parameters over a time scale of several days, which is owing to the probable apsidal 
precession or nodal precession period of the accretion disk. Suitable data sets for searching 
periodicity are the 10-day consecutive RXTE observations in 1996 and 1998. For these two 
10-day consecutive observation, all dip parameters varied significantly over a time scale of 
2-5 days (see Fig. [3] and Fig. Hj). Some of the parameters show strong periodic variations, 
especially for the dispersion in 1998 observation. The 1998 data dip dispersions exhibit 
a sinusoidal- like modulation with a period of (4.87 ± 0.14)d obtained by the sinusoidal 
fitting (see Fig. HI). This i s cons istent with the results on marginal dip phase variation 



proposed by iHomer et. al. I (120011 ). We further folded the dispersions of the entire 1998 
RXTE observation data with 4.87d period and found that all observed dip dispersions in 
1998 are coherent (see Fig. E]). It implies that the ~4.87d modulation may stably last about 
two months. A ~4d period vari ation in the dispersion, close to the ~3.9d period variation 



proposed by lChou et. al. I (120011 ) . can be seen in the 10-day consecutive observations in 1996, 



although the periodicity is not as strong as the dispersions in the 1998 data (see Fig. [3]). 

We also searched both data sets for the 3.9d and 4.87d beat side bands around the orbital 
period through the folding period search. This was provided by the FTOOLS subpackage 
efsearch. The 3.9d sidebands are clearly seen and the 4.87d sidebands are marginally detected 
in the 1996 observation. The detection of 4.87d sidebands is consistent with the negative 



superhump proposed by iRetter et. al. I (|2002[ ) . Interestingly, no clear sideband was found 
for the 1998 data set although the dip dispersion shows such periodicity. 



3.3. Orbital Ephemeris 

We collected all available dip center times from the X-ray light curves ob served from 



1978 to 2002 and applied the hnear ephemeris proposed by IChou et. al. I (120011 ): 



T^ipcenter = M J D {T D 3)50123.00944 ± 1.4 x 10"^ + 3000-6508 ± 0.0009 ^ ^ 

OD4UU 

The dip phases are first re-binned yearly and a ±0.035 phase error was assigned for phase 
jitter (see §3.41). The long-terra phase evolution is shown in Fig. [6l The linear ephemeris 



provided by IChou et. al. I (120011 ) is clearly no longer fit for 24 years of phase evolution. The 
polynomial fittings for the dip phase and time indicate that the quadratic fitting {xl = 
0.75, dof=ll) is better than the linear fitting {xl = 2.61, dof=12) (see Fig. [6]) with a 
confidence level greater than 99% by the F-test. This implies that the orbital period changed 
significantly with time from 1978 to 2002. No higher order term is required for the dip phase 
evolution because the F-test confidence level in comparison with cubic {xl = 0.72, dof=10) 
and quadratic fittings is only 74%. We therefore updated the orbital ephemeris in a quadratic 
form: 

Tn = MJD(TDE)50123.00873±0.0004+ ^°°°'^^^^^°'°°°'^ xAr+(2.67±0.56)xl0~^^xAr^ 

86400 

and the period derivative is: 

p 

^ = (1.62 ±0.34) X lO^V"^ (6) 

Porb 

This positive value means that the orbital period increases with time. §4.21 discusses further 
implications of orbital change. 



3.4. The Stability of Dips 



When the accretion stream impact the accretion disk, a bulge forms near the impact 
region. However, the stream could penetrate the a ccretion disk, a nd fo rms a ring-like en- 
hancement of the surface density (hereafter "ring"). iFrank et. al. I (119871 ) proposed a model 
of two "bulges." One lies on the impact region of the inflow stream at the accretion disk 
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edge, at about 0.8-1.0 orbital phase relative to the companion star. Another bulge, make 
up with two-phase medium of cold clouds and hot intercloud gas forms after the point of 
accretion stream and ring, is lying on the disk ring between phase 0.3 to 0.8. This model 
successfully explains the dips for the eclipsing LMXB EXO 0748-676. Most dips of EXO 
0749-676 co ncentrated around p hases 0.65 and 0.9, which correspond to the two bulges of 



this model (IFrank et. al. 1119871 ). X1916-053 has no eclipse to configure the companion star 



position as a reference point so that it is hard to verify the bulge locations responsible for 



primary and secondary dips (ISmale et. al. Ill988l ). However, we believe that the bulge near 



the edge of the accretion disk is more stable than the inner bulge because it is closer to the 
companion star, and its location in the orbital frame is less affected by the accretion disk. 

To verify the bulge locations responsible for the primary and secondary dips, we tested 
the phase fluctuations (jitters) for both kinds of dips. The dip phases folded by the best 
quadratic ephemeris obtained by primary dips show that the standard deviation of the 
primary dips (0.035, 128 dips) is significantly smaller than the secondary dips (0.073, 48 
dips) with a F-test null hypothesis probability of 6.3 x 10"^^ (see Fig. [Tj). However, since 
the quadratic ephemeris is yielded from minimizing the deviation of the primary dips (i.e. 

fitting), we also folded the dips with the best quadratic ephemeris yielded from the 
secondary dips to check consistency. The standard deviation of secondary dips (0.074) is 
still significantly larger than the primary dips (0.055), with a null hypothesis probability of 
0.005. As a result, we conclude that the primary dip is due to the budge near the edge of 
accretion disk whereas the secondary dip is caused by the budge on the ring. The standard 
deviation of primary dips (0.035) with respect to the best ephemeris gives us an estimation of 
phase jitter. This estimation is used as the systematic error of the long-term phase evolution 
in D 



4. Discussion 
4.1. Nodal Precession and Negative Superhump 



The negative superhump signal of nodal precession was detected by iRetter et. al 



(I2OO2I ) in 1996 RXTE observations. They concluded that the signal owing to persistent 
state modulation rather than the dip width (or shape) variation because the signal still ex- 
ists in the Lomb-Scargle power spectrum even when the dips were removed from the light 
curve. In our analysis, although the periodicities of dip parameters are weak in 1996 RXTE 
observation, the periodicity of dip width variation in 1998 RXTE observation is strong. 
This 4.87d period, confirmed in §3.2[ is consistent with Retter's result and provide another 
evidence of negative superhump. 
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For a bulge lies on the outer edge of accretion disk, it is believed that the bulge will 
be more opaque and larger width near the disk plane. Once while the accretion disk has 
retrograde nodal precession, the angle between di sk plane and o bserver's line of sight varies 
with nodal precession period (see Fig. 6.18 in iHellier I (I2OOII )). As a result, dip width 
variation can be interpreted easil y by the variat i on of the angle of disk plane. The 4.74d 
dip phase variation proposed by iHomer et. al. I (I2OOII ) may be induced by the dip width 
variation due to asymmetric bulge geometry. 

The mas s ratio of X19 16-053 can be interpreted by the negative superhump period. 
Chou et. al. I (1200 ll ) derived its mass ratio of 0.022 or 0.011 from its 3.9d positive super- 
hump or a Roche-lobe filled white dwarf respectively. On the other hand, for the negative 
superhump, iLarwood et. al. I (119961 ) proposed that the relation between the nodal precession 
period of the disk and the particle frequency of the outer accretion disk can be represented 
as: 

^r,. 15 , . 



Vt{R) 



15 3 . 

qr coso 

32^ 



where ujn is the frequency of nodal precession, 5 is the tilted angle of accretion disk, VL{R) is 
the particle frequency of outer accretion disk, and a ratio of specific heats 7 = 5/3 is assumed. 
For VL{R) = Suorb and assuming a smal l S, the mass ratio c an be estimated as g ~ 0.045, 
slightly larger than the ones proposed by lChou et. al. I (120011 ). All the models imply that the 
companion's mass is less than the lower limit of a normal main-sequence star (< O.OSMq), 
consistent with the previous prediction that the secondary must be a fully-degenerated or 
semi-degenerated white dwarf. 



4.2. Orbital Period Change 



Table m lists the five LMXBs with significant detection of orbital period changes prior to 
this research. Four of them (EXO 0748-676, 4U 1820-30, X1822-3 71 and Cyg X-3) have pe- 
riod derivatives inconsist ent with the standard model proposed by lRappaport et. al. I (119871 ) 



see 



Tavani et. al. II 19911 ). The standard model, in which mass loss and orbital period change 



due to gravitational radiation, predicted a positive orbital period derivative for the LMXBs 
with degenerate companions. X 1916-053 is a LMXB composed of a neutron star and a 
white dwarf. Using the model proposed by iRappaport et. al. 1 (Il987h . the orbital period 
derivative would be Port/ Port = 5.96 x 10~^^yr~^, which is a factor of 10^ to 10'^ smaller 
than the observed value. The period derivative value predicted by this model is consistent 
with the observed value only when the mass transfer is extremely not conservative (near the 
singularity in this model). 
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Tavani et. al. I (Il99ll ) proposed a radiation-driven model to explain inconsistencies be- 



tween the standard model and the observed value. They argued that the companion star 
may be illuminated by a relatively large flux of radiation from the primary star or accre- 
tion disk. Such radiation can drive a strong evaporative wind from its outer atmosphere. 
In this model, the irradiated companion can transfer mass even if it does not fully fill its 
Roche lobe. They derived an average self-sustained value of radiation-driven mass loss 
rate of — m2 ~ 10~^MQ?/r"^, which is much larger than that of gravitational radiation 
{—rriGR ~ 9.15 X IO^^'^MqUt^^). For X 1916-053, the period derivative value can be in- 
terpreted by this model with 60 — 90 % of com panion mass loss outf l ow fro m the system. 
The mass outflow has been detected in LMXBs. IChakrabarty et. al. I (120031 ) reported that 



a strong He I 1.083 /im emission line with P Cygni profile in GX 1+4/V2116 Oph through 
infrared observation. From the blue edge of this profile, they inferred that there is an out- 
flow with a velocity much faster than a typical red giant wind from this binary system. 
For X1916-053, its outflow could be verified if the P Cygni profile were detected in its op- 
tical counterpart. Unfortunately, the optical counterpart of X1916-053 is too dim (V=21) 
so that no optical spectro scopy has ever been r eported. From the mass accreted onto the 
neutron star predicted by iTavani et. al. I (jl99ll ). the X-ray luminosity can be estimated as 
1 — 4 X lO^^ero ■ s~^- This is slightly larger than the observed value (0.5 — 1.44 x IQ^^ erg ■ 



Bloser et. al. I (120001 ) ). 
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Fig. 1. — Definition of dip center times. This is an example light curve taken by RXTE in 
1996. Where /q is the persistent state fitting result, li is the observed count rate, and U is 
the corresponding time. Iq — U is the weighting factor in our definition of dip parameters. 
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Dip parameters vs the choice of boundaries 
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Fig. 2. — A test to show the dependence of dip parameters and hght curve boundaries for 
evaluation, (a) The hght curve for the test. A dip with comphcated profile appeared between 
~21300s and ~22200s (from MJD(TDB) 52542). We set the left boundary fix at point A 
(21200s, the dashed line) and the right boundary (dotted line) gradually moved from point 
B to D. Plot (b) (c) and (d) are the dip parameters calculated with different right boundary. 
It is evidentially that as long as the right boundary hes on persistent state (beyond point 
C) , the calculated parameters are insensitive to the location of right boundary. 
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Fig. 3. — The dip parameters in 1996 10-day consecutive observations. The first three panels 
represent the three dip parameters of dip with unit in seconds, while the last panel represents 
the persistent count rate near each dip. All the parameters varied over a period of several 
days, but significant periodicity cannot be detected. 
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Fig. 4. — The dip parameters in 1998 10-day observations. The most significant detection 
of periodicity is in the dispersion. Through the sinusoidal fitting (dashed curve), the period 
can be estimated as (4.87 ± 0.14) days. 
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Fig. 5. — The dispersion of the complete 1998 data folded with 4.87d period. The points 
on the diagram are 10-day consecutive observation data, whereas the plus signs are other 
observations in 1998. 
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Fig. 6.— X-ray dip phase evolution of X1916-053 from 1978 to 2002. The folding period 
is 3000.6508s and the phase zero epoch is MJD(TDB) 50123.00944. The thick curve is the 
quadratic fitting result. 



-19- 




Fig. 7. — The phase distribution of primary dips and secondary dips with respect to the 
quadratic ephemcris derived from the primary dips. The standard deviation of primary dips 
is significantly smaller than that of secondary dips. This means that the primary dips are 
more stable than secondary dips. 
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Table 1: The historical X-ray observations of X1916-053 



Year observatory observation date number of dips 

2002 XMM-Newton Sep. 25^ 

RXTE Sep. 25 7 

2001 RXTE May 25, Jun 16, 17, 30, Jul 01, Oct 01 

BeppoSAX Oct 01-02 27 

1998 RXTE Jun 24, Jul 17-26^ , Aug 01, 10, Sep 14, 16 26 

1997 BeppoSAX Apr. 27-28 9 

1996 RXTE Feb 10, Mar 13, May 05, 14-23^ , Jun 01, 

Jul. 15, Aug 16, Sep 06, Oct 29 27 

1993 ASCA May 02-03 5 

1992 ROSAT Oct 17-19 1 

1990 Ginga Sep 11-13 2 

1988 Ginga Sep. 09-12 2 

1985 EXOSAT May 24, Oct 13 10 

1983 EXOSAT Sep 17 8 

1980 Einstein Oct 11 2 

1979 Einstein Oct 22 1 

1978 OSO-8 Apr 07-14 1 



^This column represents number of dips of this year. However, some light ciirves contain no complete dips or 
have relative lower significance are folded into one fold light curve so that we treat them as one dip per fold 
light curve. 

^The simultaneous observation of RXTE and XMM-Newton 
^The 10-day consecutive observations in 1998 and 1996 
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Table 2: Dip parameters vs energy bands. The dip center time is measured from MJD(TDB) 



dip number 


Energy band (keV) 


dip center time (s) 


dispersion (s) 


equivalent width (s) 


1 




±ooiy. ozo 




AAA "Rns 




7 17 /"cirvff^ 

u./-i./ (^sorcj 


iODZ4.00i 


1 Q/l 771^ 
iy4. MO 


A'A^ /I an 




1. /-O.U \^nara j 


iODiD.OOi 


1 Q9 Q/l /I 


4oi.04 / 









187 491 


9Q1 581 




0.7-1.7 (soft) 


18766.314 


191.611 


341.541 




1.7-8.0 (hard) 


18769.707 


183.034 


250.401 


3 


0.7-8.0 (total) 


21771.504 


223.684 


415.795 




0.7-1.7 (soft) 


21776.967 


220.304 


471.319 




1.7-8.0 (hard) 


21765.254 


227.488 


370.354 


4 


0.7-8.0 (total) 


24769.473 


243.336 


474.337 




0.7-1.7 (soft) 


24771.145 


245.808 


541.436 




1.7-8.0 (hard) 


24769.922 


239.112 


417.391 


5 


0.7-8.0 (total) 


27723.248 


250.468 


457.142 




0.7-1.7 (soft) 


27729.359 


249.732 


516.507 




1.7-8.0 (hard) 


27716.568 


249.946 


404.874 



Table 3: Measured dip center time between XMM-Newton MPS detector and RXTE PGA 



dip number 


dip center time(s) 


dispersion (s) 


1 (XMM-Newton) 


15619.323 


193.48 


(RXTE) 


15603.40 


182.708 


2 (XMM-Newton) 


21771.504 


223.684 


(RXTE) 


21730.94 


225.89 


3 (XMM-Newton) 


27723.248 


250.468 


(RXTE) 


27685.86 


254.65 
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Table 4: LMXBs characteristics with measured orbital period change 



Object Name 


modulation * 


Porb (hr) 


Porbl Porb (jr ^) 


Reference 


4U 1820-30 


M 


0.19 


-3.74 X 10-s 


Chou & Grindlav r2001) 


X 1916-053 


D 


0.83 


1.62 X 10-^ 


this work 


EXO 0748-676 ^ 


D,E 


3.82 


2.7 X 10^^ 


Wolff et. al. r2002) 


Cyg X-3 2 


M 


4.82 


1.05 X 10-^ 


Sineh et. al. f2002) 


X 1822-371 


PE,D 


5.57 


3.4 X 10-^ 


Hellier et. al. ri990) 


Her X-1 


E 


40.8 


-1.32 X 10^^ 


Deeter et. al. ri991) 



*E: total eclipse, PE: partial eclipse, D: periodic dips, M: other modulation 

^The period derivative value of EXO 0748-676 came from the quadratic fitting of entire dataset from EX- 

OSAT(1985) to RXTE(2000). 
^Cyg X-3 was identified to be a High Mass X-ray binary because the companion star was confirmed as a 

Wolf-Rayet helium star with a strong wind. 



